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Abstract 

In the last years quantum correlations received large attention as key ingredient in advanced 
quantum metrology protocols, in this letter we show that they provide even larger advantages when 
considering multiple-interferometer setups. In particular we demonstrate that the use of quantum 



a. 
< 

^ correlated light beams in coupled interferometers leads to substantial advantages with respect to 



classical light, up to a noise-free scenario for the ideal lossless case. On the one hand, our results 
prompt the possibility of testing quantum gravity in experimental configurations affordable in 
current quantum optics laboratories and strongly improve the precision in "larger size experiments" 
such as the Fermilab holometer; on the other hand, they pave the way for future applications to 
high precision measurements and quantum metrology 

PACS numbers: 42.50.St, 42.25.Hz, 03.65.Ud, 04.60.-m 



The dream of building a theory unifying general relativity and quantum mechanics, the 
so called quantum gravity (QG), has been a key element in theoretical physics research for 
the last 60 years. Several attempts in this sense have been considered. However, for many 
years no testable prediction emerged from these studies, leading to the common wisdom that 
this kind of research was more properly a part of mathematics than of physics, being by 
construction unable to produce experimentally testable predictions as required by Galilean 
scientific method. In the last few years this common wisdom was challenged [IHH]. More 
recently, effects in interferometers connected to non-commutativity of position variables [71 [8] 
in different directions have been considered both for cavities with microresonators [4j and 
two coupled interferometers [5j , the so called "holometer" . In particular this last idea led to 
the planning of a double 40 m interferometer at Fermilab [9j . 

Here we consider whether the use of quantum correlated light beams in coupled inter- 
ferometers could lead to significant improvements allowing an actual simplification of the 
experimental apparatuses to probe the non-commutativity of position variables. On the 
one hand, our results demonstrate that the use of quantum correlated light can lead to 
substantial advantages in interferometric schemes also in the presence of non-unit quantum 
efficiency, up to a noise-free scenario for the ideal lossless case. This represents a big step 
forward respect to the quantum metrology schemes reported in literature [lOHlBJ , and paves 
the way for future metrology applications. On the other hand, they prompt the possibil- 
ity of testing QG in experimental configurations affordable in a traditional quantum optics 
laboratory with current technology. 

The idea at the basis of the holometer is that non-commutativity at the Planck scale {Ip = 
1.616 X 10~^^ m) of position variables in different directions leads to an additional phase noise, 
referred to as holographic noise (HN). In a single interferometer I this noise substantially 
confounds with other sources of noise, even though the most sensible gravitational wave 
interferometers are considered [5], since their HN resolution is worse than their resolution 
to gravitational- wave at low frequencies. Nonetheless, if the two equal interferometers Xi 
and X2 of the holometer have overlapping spacetime volumes, then the HN between them is 
correlated and easier to be identified [5j. Indeed, the ultimate limit for holometer sensibility, 
as for any classical-light based apparatus, is dictated by the shot noise: therefore, the 
possibility of going beyond this limit by exploiting quantum optical states is of the greatest 
interest [lU [H [E] . 
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FIG. 1: Sketch of the holometer. a) The two involved interferometers, X]^^ A: = 1,2, have input 
modes a/^ and 6/^ and output modes C]^ = Ck{(l)k) ^^d dk = dk{(j)k)^ where two detectors are placed 
for measuring the number of photons Nck(^k) ^^d Ndj^{(j)k)^ respectively, b) The interferometers 
are set in the configuration "||" and "±" according to the picture. The input modes 6/^, fc = 1,2, 
are always excited in two coherent states labeled as " |COi7)i|COi7)2", while the modes ak are 
excited in two uncorrelated squeezed vacua labeled as ^^\SQ)i\SQ)2^ [configuration (SQ)], or in a 
maximally entangled two-mode squeezed vacuum marked as "|TVFS)i2" [configuration (TWB)]. 



In the past the possibility of exceeding shot-noise limit in gravitational-wave detectors 
was suggested [161 E| ^^d, recently, realized [18] by using squeezed light. As shown in the 
following, this resource can indeed allow an improvement of holometer-like apparatuses as 
well. Nonetheless, in this case, having two coupled interferometers, the full exploitation of 
properties of quantum light, and in particular of entanglement, may lead to much larger 
improvements. 

In general, the observable measured at the output of the holometer may be described by 
an appropriate operator (7(01, 02), 0/c being the phase shift (PS) detected by X/^, /c = 1,2, 
with expectation value ((7(01, ^2)) = Tr[pi2(7(0i, 02)], where pi2 is the overall density matrix 
associated with the state of the light beams injected in Xi and X2. 

In order to observe the HN, one should compare [5j ((7(0i,02)) in two different exper- 
imental configurations of Xi and X2, namely, parallel, "||", and perpendicular, "_L" (Fig. 
[1]). In configuration "II", the interferometers are oriented so that the HN induces the same 
random fluctuation of their PSs, leading to a substantial correlation between them, since 
they occupy overlapping space-time volumes [5l[Wj. Thus, by measuring the correlation of 



the interference fringes, one can highlight the presence of the HN. Configuration "_L" serves 
as a reference measurement, namely, it corresponds to the situation where the correlation 
due to HN is absent, since their space-time volumes are not overlapping jSJ [I9j, in other 
words, it is equivalent to the estimation of the "background". The statistical properties of 
the PSs fiuctuations due to HN may be described by a suitable probability density function, 
/x(0i5 ^2), X =11, _L. In turn, the expectation of any operator 0(01, 02), or function of the 
PSs, should be averaged over fx^ namely, 

(6(01, <^2)) ^ £. [6(<^i, h)] = y (6(01, 02)) /x(01, 02) d0i d02. (1) 

As in the bolometer the HN arises as a correlation between two phases, the appropriate 
function to be estimated is their covariance in the parallel configuration ^^y [50i502], where 
50A; = 0/c — 0/c,Oi and 0A;,o are the mean PS value measured by X^, k = 1,2. Since the 
holographic noise is supposed to be small, we can expand the C operator in terms of small 
fluctuation 50^. According to Eq. ([!]) we are able to directly relate the covariance of the 
PSs to the observable quantities (see Sup. Mat. Sec. I for details): 
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(92^,^^^(01,0, 02,O)) 

Eq. ([2]) states that the covariance can be estimated by measuring the diflPerence between 
the expectation value of the operator C in the two configurations. Thus, this diflFerence 
represents the measured signal, while the coefficient at the denominator is the sensitivity. 
One has to reduce as much as possible the uncertainty associated with its measurement: 
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variances derives from the independence of the two measurements configurations. Thanks to 



the same expansions leading to Eq. (2), we can write Var^^ 
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where Var (^(^i^o, 02,o) = (C'((/>i,o, 02,o)^) — (C'((/>i,o, 02,o))^ does not depend on the PSs 
fluctuations due to the HN, but it represents the intrinsic quantum fluctuations of the 
measurement described by the operator (7(01, 02) and depends on the optical quantum states 
sent in the holometer. In particular, our aim is to look for a suitable choice of quantum 
optical state pi2 and an operator (7(01,02) that reduces this zero-order contribution to the 
uncertainty. In the following we will demonstrate that the use of quantum resources, like 
squeezing or, much more, entanglement, provides huge advantages in terms of achieved 
accuracy with respect to classical light. 

As a flrst example we consider a conflguration (SQ) where the two input modes of 
each interferometer X/^, k = 1,2, are excited in a coherent state and a squeezed vac- 
uum state with mean number of photons fik and A/e, respectively (see Fig. [l]). Since 
the difference of the number of photons in the two output ports of each interferometer, 
Nk-{(t)k) = A^cfc(0fc) — ]^df^{(/)k)^ can be used to estimate the corresponding 0/^ with sub- 
shot-noise resolution [TOl [H] [16], then reasonably the covariance ^||[50i502] can be efli- 
ciently evaluated from the covariance between A^i_(0i) and A^2-(02)- Therefore we deflne 
C(0i,02) = AiVi_(0i) AiV2_(02), with ANk-{^k) = Nk_{(t)k) - £ ^-(0^)1 (we note that 
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Nk-{4^k) = ^± Nk_{(j)k) = £ Nk-{4^k) ) as consequence of the properties of /a:(0i, 02), 



see Sup. Mat. Sec. I). 

Fig. [2] shows the corresponding uncertainty at the zero-order given in Eq. ^\ assuming 
identical input states [jik = M ^^d Xk = X^ k = 1,2), the minimum is achieved for 0i^o = 
02,0 = 7r/2, and reads (see Sup. Mat. Sec. II): 



In perfect analogy with the PS measurement for a single interferometer p!6l [T7] . if /x ^ 
A :» 1, then we have the optimal accuracy U^q ^ (2^/2A/x)~-^. This represents a strong 
advantage in terms of uncertainty reduction [of the order (4A)~^] with respect to classical 
case UqI ^ a/2/m, i-e. when only coherent states are employed. Nevertheless, an impor- 
tant difference between the single interferometer PS measurement, involving a flrst order 
moment of the photon number distribution, and the covariance estimation, involving the 
second order momenta, arises: while in the flrst case the uncertainty scales as the usual 
standard quantum limit one, oc /x~^/^, in the second case it scales oc /x~^ (neglecting the 
little relative contribution of the squeezing to the intensity) . We remark that the advantage 
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FIG. 2: The uncertainty of the phase covariance log^^o^s-o when the holometer is fed by two 
independent squeezed states plus two coherent fields. Here ^i^o — 02, o — 00 is the central PS of 
the interferometers, Ai = A2 = A is the intensity of each squeezing, /xi = /X2 = /i is the intensity of 
each coherent beam, and their phase difference is set to zero, (a) /x ^ A: when the coherent and 
squeezed beams have similar intensities the noise reduction is lower bounded, (b) /x » A: in this 
regime quantum noise reduction increases with the amount of squeezing, offering a strong noise 
suppression if high level of squeezing can be reached. The region of the minimum runs at 00 = 7r/2. 



of the present scheme is based on the independent improvement of the resolution of each 
interferometer which is itself limited by the amount of squeezing (see Sup. Mat. Sec. II). 
However, the aim of the holometer is to couple Xi and X2 minimizing the noise on their out- 
puts correlation, namely regardless of the noise in the single interferometer. This suggests 
that quantum correlated states, coupling Xi and X2, could further enhance the performance 
of the holometer. 

To this aim, we consider a new configuration (TWB) where modes ai and a2 of Figjl] 
are excited in a continuous variable maximally entangled state, i.e. a two-mode squeezed 
vacuum state or twin-beam state, \TWB))a^^a2 = 5'i2(C)|0)ai,a2 = Y.n=o^n{C)\n)ai\n)a2^ 
where 5^12 (C) = exp(C a\al — C* aia2) is the two-mode squeezing operator. This state can be 
easily produced experimentally, for example by the parametric-down conversion process [29j . 
If we set ( = |C|e^^^ and introduce the mean photon number per mode A = sinh^ l^, then 
CniO = (l + A)-i/2 ^^i ^ ^-1^ g-z2^c]"^/^ [22J. The input modes 61 and 62 are still excited in 



two coherent states, so that the four-mode input state is 
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One of the peculiarity of the state |TWB))ai,a2 is the presence of the same number of 
photons in the two modes p!4l [23l - [25] . then each power of the photon number difference of 
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the two modes is identically null, «,,«2((TWB|(iV«, - Na,)^\TWB))a,,a2 = 0, VM > 0. 

We also observe that, in the absence of the HN and choosing the optimal working regime 
0/c,o = 0, A; = 1,2, Xi and X2 behave like two completely transparent media for their input 
fields [see again Sup. Mat. Sec. II, Eq.s (4)]. In particular output modes Ci and C2 exhibit 
perfect correlation between the number of photons, which directly comes from the input 
modes ai and a2, leading to the the natural choice of the observable as the fluctuation of the 



photon numbers difference, (7(01, 02) = A^ 



iVci - N,, 



Indeed, the numerator of Eq. (4), 



Var < A^ M^ci(O) — Nc^iO) > is identically null, while the denominator reads: 






Nc^iM - Nc,{<t>2,o)] I ^) = -^VA(l + A)/icos[2(^c " ^a)] (6) 



that is non-zero for both A, /x 7^ and it is maximized for 9(^ — 9a = 0. This quantity 
represents also the coefficient of proportionality in Eq. ^ between the covariance of the 
HN and the measured signal. It is worth noting that, even though for (1)^,0 = the coherent 
state gives no contribution to the output modes Ci and C2, being completely transmitted 
to the complementary modes di and ^2, when ffuctuations of the PS occurs a little portion 
will be reflected to the monitored ports and this guarantees the sensitivity to the HN PSs 
covariance. 

Thus, the correlation property of the TWB state leads to the amazing result that the 
contribution to the uncertainty coming from the photon number noise shown in Eq. Q is 
^TWB ^ (when A, /x 7^ 0), representing an ideal accuracy of the interferometric scheme to 
the PSs covariance due to HN and the main achievement of the present study. 

The question that now arises is how and at which extent our conclusions are affected by 
a non unity overall transmission-detection efficiency r] (see Sup. Mat. Sec. HI). In Fig. [3] we 
plot Z//(^^ for the SQ, TWB and classical coherent state (CL) approaches, as a function of 
77 (assumed to be the same for both the interferometers) for a modest level of non-classical 
resources (A = 0.5). As one may expect, SQ exhibits a little advantage in this regime. 
However, in the high efficiency region (albeit with values reasonable achievable with current 
technologies) the TWB-based approach provides a signiflcant improvement not only with 
respect to classical set-up, but also respect to SQ. 

Focusing on the limit of very small quantum resources, i.e. A <C 1 and /x 3> 1, then 
^sq/^cl ^ 1 ~ 2r]\/\ and Z//twb/^cl ^ \/2(l — r])/r]. For a small amount of squeezing 
the quantum noise ZYgg not surprisingly approaches the classical case, while for TWB, we 
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have a degradation of the performances with respect to the ideal case (rj = 1). Anyway, an 
improvement with respect to the classical case is kept for rj > 2/3, demonstrating that a 
relatively faint TWB can provide an interesting improvement in the HN detection. 

In the opposite limit of high quantum resources exploited, i.e. /x ;» A :» 1, Z^/sq/^cl ^ 
(1 — 7y) +7y/(4A) and Z//twb/^cl ^ 2a/5 (1 — rj) reveal that the performances of the quantum 
strategies are limited by the presence of the terms (1 — ry). Here the main difference between 
SQ and TWB is that for rj ^ 1 SQ exhibits an uncertainty lower bounded to (4A)~^, i.e. 
depending on the squeezing intensity. On the other hand, TWB approach beats the classical 
one for rj > 0.683, while foirj^l it goes to zero: this demonstrates that the use of quantum 
light can largely improve the performances of a holometer addressed to test quantum gravity 
models. 

A last source of noise, which could affect our results, may derive from the radiation 
pressure (RP) [lEl Ell [28] . However, our model is perturbative in phase fluctuations 50/^ and, 
according to Eq. ([5]), the main noise contribution should come from the zeroth-order term 
corresponding to the photon noise. RP is assumed to introduce a second-order contribution 
to the uncertainty [see Sup. Mat. Sec. IV for details], that is related to the light fluctuation 
in the arms of the interferometers and to the phase shift induced by the mirrors recoil (the 
latter is given by {ujT/2mc)V^ where u) and V are the central frequency of the light and the 
momentum of the photon, respectively, r the measurement time and m the mirrors mass 
[16]). In the case of a single interferometer fed by squeezed light, the amplitude of the RP 
noise varies with the squeezing parameter at the opposite of the photon noise, namely, if 
one decreases the other increases and viceversa, thus an optimum regime must be found. In 
the context of our proposal the behavior is similar, but for reasonable values of the involved 
parameters RP noise is completely negligible (justifying our perturbative approach) . Fig.|3] 
shows how RP noise starts to affect the uncertainty for an average coherent field intensity 
M ^ 10^^ for T = 10~^ s (the photons introduced by the quantum modes are negligible), 
that is a power larger than hwji/r = 3.3 x 10^ Watts! Since the HN must be seek in the 
region of the MHz, i.e. for a short measurement time (r ^ 10~^ s), the RP contribution 
would be significant for power values larger than 10^^ Watts, a value well far away from 
the current and future interferometry technology. One can be surprised that in the present 
scheme the radiation pressure is negligible, while it is not always the case for the standard 
phase measurement involving a single interferometer. However, we stress that here we are 



measuring a phase covariance between two interferometers, instead of the phase values in 
single ones. 
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FIG. 3: Uncertainty reduction normalized to the classical limit U^j^ [26] . The solid lines represent 
the uncertainties only due to the photon noise, corresponding to the zero-order contribution [see 
Eq. Q]. The coherent field intensity is set to /x = 2 x 10^^, while the twin beams and independent 
squeezed beams intensities are Ai = A2 = A = 0.5. The dashed lines represent the second- 
order uncertainties including the RP contribution. For its calculation the measurement time is 
chosen r = 10~^ s, the mirror mass m — 10^ Kg, and the central angular frequency of the light 
uj = 3.14 X 10^^ Hz (a wavelength of 600 nm.) 



In conclusion, besides our analysis concerning the use of two independent squeezed states, 
which substantially confirms the advantages already demonstrated for gravity waves detec- 
tion, the investigation carried out exploiting entanglement leads to the unprecedented result 
that, in an ideal situation, the background noise can be completely washed out. This 
achievement not only paves the way for reaching much higher sensibility in the bolometer 
in construction at Fermilab or for the realisation of a table top experiment to test quantum 
gravity, but also sheds some first light on new unexpected opportunities offered by the use 
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quantum states of light for a fundamental reduction of noise in interferometric schemes. 



I. SUPPLEMENTARY MATERIAL 



A. Derivation of Eq. (2) of the main text 

Since the fluctuations due to the holographic noise (HN) are expected to be extremely 
small, we can expand (7(01, (^2) around the phase shift (PS) central values ^i^o, 02,0, namely: 

C{(t)l, h) = C^(01,O, 02,o) + ^k 5<^,C(01,O, 02,o) S(f)k + -S/e 9j^,^^C(01,O, 02,o) S(f)l 

+a5,,^,C((/)i,o, 02,o) ¥1^2 + 0(¥') (7) 

H — 0^,0, and d^h^kC{(j)i^o^ (/)2,o) is the {h + k)-th. order derivative of 0(01, 02) 



^M' ^'^ 



1,2 



where 50/^ 
calculated at (pk 

In order to reveal the HN, the holometer exploits two different conflgurations: the one, 
''II", where HN correlates the interferometers, the other, "± ", where the effect of HN 
vanishes (Fig. (1) of the main text). The statistical properties of the PS fluctuations 
due to the HN may be described by the joint probability density functions /y (01,02) and 
/±(0i5 02)- We make two reasonable hypotheses about /x(0i,02), x =||,-L. First, the 
marginals J'x {(j^k) = /d(/)/j/a;((/)i, ^2), h^k = 1,2 with h ^ k^ are exactly the same in 



rW/ 



rW/ 



can write /±((/)i, ^2) = J='^^\(t)i)J^^^^ 



the two conflgurations, i.e. J'n'^\(j)k) = J~'±\4^k)' one cannot distinguish between the two 
conflgurations just by addressing one interferometer. Second, only in conflguration "±" we 

2), i.e., there is no correlation between the PSs due to 
the HN. Now, the expectation of any operator 0(01, ^2) should be averaged over fx^ namely, 
(0(01,02)) -^ Sx 0(01,02) = /(O(0i,02)) /x(0i,02) #1 #2- In tum, by averaging the 
expectation of Eq. ([7|, we have: 



^cr. 






(8) 



where we used £x[5(l)k] = 0. Then, according to the assumption on fx{4'i,<f>2) we have 
^ll [S(j)l] = £^ [5(f)l] and £^ [(501(502] = E^ [8(t)i] £± [S(j)2] = 0, and from Eq. g foUows that 
the PS covariance may be written as in Eq. (2) of the main text, namely: 
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that is proportional to the difference between the mean values of the operator (7(01, 02) as 
measured in the two configurations "|| and "_L". 

Indeed, one has to reduce as much as possible the uncertainty associated with its mea- 
surement, which reads as: 



U{5(j)i5(j)2) 
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Under the same hypotheses used for deriving Eq. ([9]) we can calculate the variance of 
C(0i,(/)2) as 



Var^ 



where: 



C(01,02 



Var 



C*(01,O, 02,O) + Sfc Afcfc E^ [502] + Ai2 e. [501^02] + c(<50') (h) 



Mk = (C(0l,o,02,o)55^,^j,C'((/)i,o,02,o)) 

+ (KC(0i,o,02,o)]') - (C(0i,o,02,o))(9|„^,C(0i,o,02,o)) 

Al2 = 2(5(01,0, 02,o)a|^,^,^(01,O,02,o)) 

+ 2(a^iC?(0i,o,02,o)<9,/,2C(0i,o,02,o)) - (C?(01,O,02,o))(5|j,,/,2C*(01,O,02,o)) 



(12) 



Analyzing expression (11), we note the presence of a zeroth-order contribution that does not 



depend on the PSs intrinsic fiuctuations, and represents the quantum photon noise of the 
measurement described by the operator (7(01,02) evaluated on the optical quantum states 
sent into the holometer. The statistical characteristics of the phase noise enter as second- 



order contributions in Eq. (11) from each interferometer plus a contribution coming from 



phase correlation between them. 

This work addresses specifically the problem of reducing the photon noise below the shot 
noise in the measurement of the HN, therefore in the following and in the main text starting 
from Eq. (3), we assume the zero-order contribution being the dominant one. Of course, 
this means to look for the HN in a region of the noise spectrum that is shot-noise limited. 
Since the HN is expected up to frequencies of tens MHz, it follows that all the sources of 
mechanical vibration noise are suppressed. The only additional source of noise could be 
the radiation pressure, which is itself related to the light fiuct nation in the arms of the 
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FIG. 4: The action of an optical interferometer, e.g., a Michelson interferometer (left), with input 
modes a and h and output modes c and d measuring an overall (j) phase shift between the arms can 
be summarized as a beam-splitter like system with transmittance cos^((/)/2) (right). 



interferometers. Nevertheless, we will demonstrate in Sec. |ID| that it is indeed negligible in 
our regime of interest. 



B. Elements for deriving Eq. (5) of the main text 

At first, we briefly review a basic aspect of the quantum description of optical interfer- 
ometers. We consider two input modes described by the bosonic fleld operators a and h 
{[a^a^] = 1 and [6,6^] = 1), entering the two input ports of the optical interferometer, for 
instance a Michelson interferometer, as depicted in Fig. [4| whose main component is a 50:50 
beam splitter (BS). The two modes interfere a flrst time at the BS, then they are reflected 
by two mirrors and interfere a second time at the BS after having gained an overall (j) phase 
shift) induced by the difference between the optical paths of the two arms. We refer to the 
outgoing modes as c and rf, respectively (Fig.ji]). The input/output relations can be written 

as [29]: 

c = c((/)) = cos((/)/2) a - i sm{(f)/2) 6, 

(loj 
d = d{(j)) = - cos(0/2) b + i sin((/)/2) a. 

We note that such input/output relations are equivalent to those of a BS with overall 
transmittance cos^ (0/2) (Fig.ji]). 

It is well known, in quantum optics community, that coherent light itself is not the 
optimal solution for PS estimation of an interferometer, and that the use of squeezed light 
may enhance the performance up to the Heisenberg limit [30]. In this section we assume 
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that the input modes ak and bk are excited in a squeezed vacuum state \Ck)ak = Saj^{^k)\0)ak 
and the coherent state \ak)bk = Afc(<^/c)|0)6fc, respectively, where Saj^i^k) = ^w[Ck (4)^ ~ 
^k i^kY] i^ the squeezing operator and D})^{ak) = exp{ak bl — al bk) is the displacement 
operator. From now on, for the sake of simplicity, we drop the subscript fc = 1,2. If we 
set ^ = \^\e^^^ and a = y//ie^^^, then A = sinh^ \^\ and fi represent the average number of 
photons of the squeezed vacuum and of the coherent state, respectively. In the case of a 
single interferometer, when PS is estimated from the measurement of the relative number 
of photons iV-(0) = Nc{(j)) — Nfi{(j)) [with Nc{^) = c^c and Nd{^) = d)d) note that c and 



d depends on 0, see Eq.s (13)], one has (A^_((/))) = {fi — X) cose/), and in the limit /x :» A, 
around the optimal working point cj) = 7r/2, the predicted uncertainty is e^'^'/v^' ^•^' below 
the shot noise or standard quantum limit [31j. 

Following the same line of thought we investigate the possibility of exploiting the squeezed 
states of light also in the case of the estimation of the PSs covariance between the two inter- 
ferometers. Since Nk-{4^k)^ A; = 1, 2, can be used to estimate the PS (/)/e of the interferometer 
X/e, it is reasonable to evaluate the covariance of the PSs S\\ [501 502] in Eq. ^ from the co- 
variance between iVi_((/)i) and iV2-(02)- We define (7(01,02) = ^Ni_{(/)k) ^^2-{(t)k)-> with 



AiV,_(0,) = Nk 



- £ 



m. 



and £ 



m. 



£\\ 



m. 



£i 



Nk-{ct^k) 



according to the properties of the marginal distribution of /x(0i,02)- Thus, the expected 
value of (7(01,02) in each configuration is effectively the covariance between Ni_{(j)k) and 



iVi_(0fc) iV2_(<^fc) 



£ 



Ni-{<l>k) 



£ 



N. 



.) 



A^2-(0fc), i-e., £x C{(j)i,4>2) = £x 

We inject the state |t/^) = |^i)ai ® |tti)6i ® 1^2)02 ® \oi2)b2 through Xi and X2 considering 
symmetry properties of the states, Cfc = ^) ct/c = ct and of the interferometers 0i^o = ^2,0 = ^0, 
and setting 6*0.^1 — 6^^i = 9a,2 — ^5,2 = 0, k = 1,2. 

According to Eq. (4) of the main text, and by substituting the input-output relations 
(13) into the definition of the observable C{(f>i,4)2) = ANi_{(f>i) AN2-{4>2) the minimum 



uncertainty becomes: 



mw ,. r- A + u(l + 2A-2VA + A2) 
W(°)(/i,A,(/>o = 7r/2) = 72 ^^ (A-/1/ ■ 



(14) 



As a final comment, we remark that the advantage of the present scheme is naturally 
limited by the amount of squeezing. In fact, by setting = 7r/2 in Eq.s (13), we observe 
that the interferometer behaves as a 50:50 BS, and the measurement of Nk- corresponds to 
the quadrature measurement of the input mode ak by the well known homodyne detection 
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technique, where the mode bk plays the role of an intense local oscillator (/x :» 1). Increasing 
the intensity A of the squeezed field in the modes ak means, by definition, reducing the 
quantum fiuctuation on the measured quadrature and thus enhancing the accuracy of PS 
estimation in each interferometer. 

C. The effect of losses 

The overall effect of losses can be modeled by means of BSs with a suitable transmittance. 



Formally, this corresponds to the substitution of the output modes c and d of Eq.s (13) with: 

C^ Crj = 77^/^ C + (1 — 77)^/^ Vr 



^ci 



(15) 
d^drj = r]^/^ rf+ (1-77)1/2 y^, 

that is the transmitted outputs of two identical BSs, both with transmittance 77, where the 
modes Vc and v^ are taken in the vacuum state |0^J ® |0^d)- Thereafter, the calculation 
of Z//gQ and ^xwb ^^ ^ function of r] is performed in complete analogy of the lossless case 



see Sec. IB for independent single mode squeezed states, while for TWB see the specific 



paragraph of the main text). 



D. Radiation pressure noise 



Our analysis is based on the reasonable hypothesis of small phase-shift fiuct nations. 
Therefore the main contribution to the noise should come from the intrinsic photon noise in 



the measurement, given by the zero-order term in Eq. (11). However, it is important to test 
this assumption including, at least, the effect of the other unavoidable contribution to the 
noise due to the radiation pressure (RP). In the context of our proposal, we demonstrate 
that, for reasonable values of the involved parameters, RP contribution is negligible. 

The statistical properties of the RP noise are determined only by the fiuctuation of 
the photon number inside the arms of the interferometers, that is independent of the 
HN. In particular, the fiuctuation can be written as the sum of the two independent 
contributions, namely, 6(l)k = 54^k,HN + 54^k,RP) while the global probability density is 
the product of the HN density function with the RP one F{(J)i^hN)4^i,hN)4^i,rp4^2,rp) = 
/x(0i,ifA^5 4^2,hn) 5'(0i,i?p, 02,i?p), where we assume that g is the same for x =||, -L. Since the 
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variance and covariance of the PSs over F are: 



^x 



I] 



^x 
^x 



k,HN\ 



Sj, [S(l)iS(l)2] = £x [Hl,HN5(t)2,HN] + £ 

respectively, the equation ^ can be written as: 



2 1 
k,RP\ 1 



£\\ 



£\\ 



£i 



h,RP0(p2,RP\ , 



Ci4>l,4>2) 



(16) 
(17) 



n,HN^^2,HN\ 



(18) 



(^|„^,C'(0i,o>2,o)) 

If the HN is absent or by far smaUer than the other sources of noise, as expected, the 
measurement uncertainty stemming from the photon noise and the RP noise can be obtained 



by using Eq. (11) in the numerator of Eq. (10): 



U^^\5ct)i5ct)2) 



\ 



Var 



C{(i)lfi, (f)2,o) + Sfc Akk £ [^(t^%Rp] + ^12 £ [S(I)1,RpH2,Rp] 



(9|,^,C(01,o,02,o))l' 



(19) 



The difference in the transferred momentum to the two mirrors of the interferometer I^ 
is Vk = {2hw/c)n-^k-> where ri-^k represents the photon numbers difference in the two arms 
and (jj is the central frequency of the light. The arms length difference induced by the RP 
can be written as Zk^Rp = {T/2m)Vk^ t being the measurement time and m the masses of 
the mirrors, and the corresponding PS is (j)k^Rp = {^/c)zk^Rp- Therefore, the fluctuation 
5(j)k^RP is proportional to 5n_,/e = ^-,/c — (^-,/c) ^^d the variance and the covariance are: 



£ 



£ [HIrp] 



1,RP 0(P2,RP\ ~ 



c^m 



huj'^T^ 



c^m 






(20a) 
(20b) 



respectively. The quantum expectation values in the right-hand sides of Eq.s (20) can be 



evaluated easily by writing n^^k = <^I out^k,out — &I out^k,out according to the beam splitter 
input-output relations ak^out = 2~-^/^(ajt + bk) and bk^out = 2~-^/^(6/e — ^k)- In the case of 
independent injected squeezed states \ip) = \^)ai ® \(^)bi ® \0a2 ® \^)b2 we have: 

{Sri^-,k)sQ = A + /x |l + 2A + 2VA(l + A)cos[2(e,,, - %,)]} (21a) 

{Sn_^iSn_^2)sQ = 0. (21b) 

where we have set the optimal working point cpi^o = 02,0 = 0o = 7r/2. We note that here 
the contribution of the covariance of the photon number difference is null, because of the 
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independence of the squeezed states sent in the two interferometers, while the variance in 
each interferometer agrees with the one reported in [31j. As in the usual treatment of a 
single interferometer fed by squeezed light, the amplitude of the RP noise varies with the 
squeezing parameter at the opposite of the photon noise (if the photon noise reduces the RP 
noise increases or viceversa) [31j. 



When the expectation values in Eq.s (20) are evaluated for the TWB case, i.e. for the 



state I?/;) = |TWB))a^^a2 ® \^)bi ® \^)b2^ we obtain: 

{Sri^_^k)TWB = A + //(1 + 2A) (22a) 



(5n_,i 5n_^2)TWB = 2^l^/X{l + X)cos[2{9^^k - ^c,^)]- (22b) 



where we have set the optimal working point ^i^o = 02,o = 0o = 0. Comparing Eq.s (21) 



and Eq.s (22) one notes that the single interferometer contribution is of the same order for 
both "SQ" and "TWB" configurations, while the covariance of the "TWB" is non null and 
basically similar to the variance of the RP of the single interferometer (if 9^^^ — ^c.k = 0, as 
considered through all the paper) . Even if the effect of the RP noise on the final measurement 



is given by Eq. (19), and requires the non trivial calculation of the coefficients Akj^ this 
suggests that the RP could be more effective in the case of TWB. We do not report here the 
lengthy calculation which leads to a cumbersome result; nevertheless, we show in Fig. [5] the 
behavior of the uncertainties U as functions of the mean number of photons of the coherent 
state ji. It is worth noting that RP contribution is negligible for reasonable value of /x, while 
for non-realistic higher value of /x where RP becomes relevant, TWB is a bit more penalized 
comparing with the case of independent squeezed states. 

Summarizing, the advantage of using twin beams is not affected substantially for reason- 
able values of the parameters, as shown also in Fig. 3 of the main text, where the uncertainty 



U^^\5(t)i5(t)2) of Eq. (19) is plotted as a function of the detection efficiency (see Sec. IC). 
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FIG. 5: Uncertainty reduction normalized to the classical limit ZY^^, versus the normalized mean 
number of photons of the coherent fields ji/ R. R — ( ^^ j is the characteristic adimensional 
parameter connected with the phase-fluctuation response to the photon number fluctuation, ac- 
cording to Eq. (20): one has i? ^ 8.6 x 10^^ for the values of r, m and uo chosen in Fig. 3 of the 



main text. The solid lines represent the uncertainties only due to the photon noise, corresponding 



to the zero-order contribution [see Eq. (19)], while the dashed lines represent the second-order 
uncertainties including the RP noise. The twin beams and independent squeezed beams intensities 
are Ai = A2 = A = 0.5 and the overall transmission-detection efficiency 77 = 0.98 
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